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; Abstract 

We model the NN system at rest as a baryon number zero lump in the Skyrme 
model including the to field. We integrate the classical equations of motion 
from the highly non-perturbative annihilation region into the non-interacting 
| radiation zone. Subsequent coherent state quantization of the radiation field 

gives a good description of the pion spectrum from annihilation at rest. 
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I. INTRODUCTION 



We have recently shown that low energy nucleon-antinucleon annihilation directly into 
pions can be well described as occurring through a coherent classical pion wave that emerges 
very quickly from the annihilation region IQ , @] • The physical, quantum pions detected from 
annihilation can be described by quantizing this classical pion burst using the method of 
coherent states ||. We project these coherent states onto states of fixed four-momentum 0] 
and definite isospin ||. A correct picture of the pion spectrum, pion number probability and 
charge averages emerges from this simple picture. Our approach is inspired by calculations 
of Skyrmion anti-Skyrmion annihilation, where it is found that the classical Skyrmion fields 
annihilate into a burst of pion radiation nearly as soon as the Skyrmion and anti-Skyrmion 
touch, and that that pion wave takes away the baryon number and energy as fast as causality 
will permit @, @. This picture, and our first treatment dealt only with the direct pions and 
did not address the roughly 40% of annihilations that go first through meson resonances, 
mostly omega and rho mesons, and then to pions [§]. The purpose of this paper is to begin 
to address this deficiency by including omega mesons in our description of annihilation. We 
will return to rho mesons in a subsequent paper. 

In our previous work ||, we approached the classical pion wave from annihilation 
phenomenologically, and concentrated on projections of the coherent state. We found that 
a very simple parameterization of the wave source gave a good description of the pion 
spectrum. The only parameter we needed was the overall amplitude of the wave, and this 
we could fit to the total energy released. When coupling to omega mesons the situation is 
more complicated, since both the relative amplitude and shape of the omega and pion waves, 
even classically, should be dynamically determined. Thus we need to make a dynamical 
theory of the coupled pion and omega fields. Again we do that classically. We construct a 
classical picture of the annihilation proceeding into interacting pion and omega fields. We 
propagate those fields, classically, into the radiation zone, where they no longer interact, and 
then quantize the radiation into coherent states for the pions and omegas separately. Finally 
we allow the omega mesons to decay into pions, assuming that that radiation is incoherent 
with respect to the direct pions. This approach of treating the non-perturbative aspects of 
QCD classically and then quantizing the subsequent radiation in the radiation zone using 
the methods of coherent states has recently appeared in a number of applications ||, JTD| . 

For a dynamical theory of coupled, classical pion and omega fields, we use the Skyrme 



model modified to include omegas ||11|| . Recall that the usual Skyrme lagrangian |12[ de- 
scribes a classical pion field only, and that baryons emerge as topologically stable configu- 
rations of that pion field. Also recall that this theory is a candidate for what QCD looks 
like in the classical or large number of colors limit. To include omega mesons in this for- 
malism, the classical omega field is coupled to the baryon current. This helps to stabilize 
the Skyrmion. To model annihilation at rest in this formalism, we follow our earlier work 
|7j and begin with a spherically symmetric "blob" of Skyrmionic matter of size about 1 fm, 
baryon number zero and the total energy of two nucleons. This initial configuration is not 
a static solution of the Skyrme equations, but rather evolves according to those equations 
into radiating and interacting pion and omega fields. As these fields radiate outwards, they 
diminish like 1/r, decouple and become radiating free fields. We use these radiation fields to 
construct coherent states for the omegas and pions. We then project these coherent states 
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onto states of good four-momentum and isospin, as before, and examine the consequences 
for nucleon-antinucleon phenomenology. We find very good agreement with the principal 
features of annihilation, now including the omega channel. The coherent state picture seems 
to work in spite of the fact that the mean number of omegas is of order one. We assume 
that omega decay occurs outside the interaction region of the classical fields and therefore 
we do not treat the pions from omega decay as coherent with the direct pions. But we do, 
of course, treat the pions from the omega as part of the observed pion spectrum. 

In Section II we present the classical Skyrme lagrangian with coupling to the omega 
field and the dynamical equations for propagation of the fields from a spherically symmetric 
source. We also show the classical field configurations that emerge from the decay of an 
initial "blob" meant to model annihilation. Section III discusses the quantization of the 
classical it and u fields using the method of coherent states modified to include projection 
onto fixed isospin and four-mometum. Section IV gives the results of this formalism for the 
pion spectrum from annihilation. We find that the mean number of pions, the variance, 



and the entire shape of the pion number spectrum closely resemble experiment ||13|| . In 
our previous, pions only, calculation [Q, we found that only an even number of pions could 
come from even isospin states and only an odd number from odd. The pions from omega 
decay remove this artifact, giving smooth pion number probabilities in each isospin channel. 
Section V presents some conclusions and directions for future work. 



II. CLASSICAL CALCULATION 

We begin with a classical theory of the coupled pion and omega fields based on the 
Skyrme model. This coupled theory was written down by Adkins and Nappi JTT[] . In terms 



of the usual SU(2) valued unitary field U of the Skymre picture the lagrangian is given by 

C = C n + Cuj + £i n t (1) 

Ar = ~ f -fTr (d,Ud^) + \mlf$Tr{U - 1) (2) 



where 



i 

U = exp(— t ■ 7?) (3) 



A, = --u^ + —u^ (4) 



C mt = (3u„B» (5) 



B» = -^e^ al3 Tr \(U j d u U)(U^d a U)(U^d p U)] (6) 
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The first term, C W} is the standard non-linear sigma model term with a pion mass term 
included. It has two paramenters, /„., the pion decay constant, which it is customary to vary 
a little to fit the nucleon mass, and the pion mass, which we do not vary. Note that there is 
no Skyrme term in the lagrangian, since the coupling to the omega now stabilizes the theory. 
The second term, is the free omega lagrangian, where ui^ is the omega field, ui^ the 
corresponding antisymmetric field strength tensor, and m the omega mass. The interaction 
term, Ci nt , couples the omega field to the baryon current, _B M , and (3 is a coupling parameter. 
We take (3 = 15.6 and f n = 62.0MeV from the work of Adkins and Nappi p| . C int describes 
the coupling of the omega to three pions. 

To model annihilation, we will take for our initial configuration a spherically symmetric 
"blob" of matter with zero baryon number but total energy of two nucleons. Thus we only 
need solutions of the lagrangian that have spherical symmetry. This is a great simplification. 
For the spherically symmetric case we may write 



U — exp(ir • rF(r,t)) 
in which case the lagrangian becomes, 



(7) 
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'dFV _ 2 sin 2 F 




+ m 2 / 2 (cosF-l) 



+ 
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(8) 



As we can see, ug and are completely decoupled from the ti field. We therefore shall not 
consider these two components. 

The energy of the system is given by 



H = 47r r 2 drH(r) 
Jo 



(9) 



where 



H(r) = f 




dr 



'dF\ 2 sin 2 F 
dr j r 2 
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+ /> 2 (l-cosF) + - 
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(10) 



The equations of motion are given by 
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dt 2 



dr 2 



2d_ 

r dr 



+ m z 



(3 sin 2 FdF 



2vr 2 



dt 



We model the initial static field configuration as 



F(r,t = 0) = h 



r 2 + a- 



exp(— r/a) 



(13) 



(14) 



with a = l/m n and h determined by the total energy. This form corresponds to a compact 
initial "blob" of zero baryon number. Since we want to begin with the system at rest we 
take 



lj = 0. 



We may also take 



(r, t = 0) = 



since B r = at t = 0. These conditions determine u>o(r,t = 0), we find 



Uo (r,t = 0) = r dr'G(r,r) \-(5r' 2 B Q {r) 
Jo L 



where 



G(r, r ) 



1 



Irarr 1 



^—m\r—r \ ^— m(r+r ) 



(15) 



(16) 



(17) 



(18) 



Using the initial values we integrate the equations of motion to determine the fields at later 
times. As the fields propagate outwards, they diminish in size, making the non-linear and 
coupling terms in the lagrangian less important. Thus we can define a radiation zone where 
the fields propagate as linear free fields. In the radiation zone, we can calculate the pion, 
f(k), and omega, g(k), momentum distribution amplitudes from the expressions 



and 



dN n 

~¥k 



d 3 k 



\g(k)t 



1 



ixkl 



J 7i 



drr 2 3l (kr)(k^ + t-^)F(r,t) 



f°° 

drr 2 ji(fcr)(fcg' + i—)u r (r,t) 



(19) 



(20) 



respectively. 

The results of our calculations are shown in the figures. We have fixed h in our initial 
configuration by the requirement that our initial state have total energy of 2 GeV. Figure 
1 shows the pion field configuration, F, as a function of r and t, both in Fermi. We see 
that it emerges as a coherent pulse. Much of that pulse travels nearly along the light cone, 
because of the light pion mass. The uo and uj r fields are shown in Figures 2 and 3. Recall 
that u) r is initially zero. We see that both emerge more slowly and with more complex 
behavior than the pion field, due to the large omega mass. The following figures show the 
energy density. First we see (in Figure 4) the total energy density then in Figures 5 and 
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6 the energy density in the pion field and in the r-component of the omega field. That is 
the part of the omega that corresponds to radiation. In all these plots we show the energy 
density multiplied by Anr 2 , that is the total energy at any time is obtained from them by 
integration over dr only. For the total, this integral must yield 2 GeV at every time, and we 
use this test to verify the stability of our calcuation. Without the factor of r 2 , the energy 
density would be seen to decrease rapidly with r, reflecting the corresponding decrease in 
the fields. This descrease is most dramatically seen in Figure 7 where we plot the interaction 
energy density. We see that past the first few time steps, the interaction energy disappears, 
signifying the decoupling of the fields. Beyond this point not only do the fields decouple, 
but, to a very good approximation, they obey a linear wave equation. This region is the 
radiation zone, and there we can use Eqs. (HD and © to obtain the pion and omega 
momentum distributions. In Figures 8 and 9 we plot the reduced momentum distribution 
densities which are defined by 

p f (k)=4nmk 2 \f(k)\ 2 (21) 

and 

p g {k) = A7rmk 2 \g{k)\ 2 (22) 

as functions of scaled dimensionless wave number k = k/m. We see that the pion distribution 
is quite sharply peaked, roughly at a wave number corresponding to the size of the initial 
distribution. The omega distribution is much broader, largely because that field propagates 
with a mass that is large compared with other scales in the problem. The momentum 
amplitudes, f(k) and g{k) are the basis for our coherent state quantization. 

It is amusing to compare the dynamical calculation reported above with our previous 
phenomenological description of annihilation into pions only, [|TJ, ||. We can repeat the 
dynamical calculation done above using the same initial form, F(r,t = 0), for the baryonic 
configuration, but this time having no omega field and using the standard Skyrme lagrangian 
T2\. This is the appropriate dynamics for our previous, pions only, picture. Once again 



we propagate the pion field into the radiation zone and calcuate the momentum density, 
f(k). The corresponding reduced distribution densities are shown in Figure 10, where it 
is compared with the phenomenological form used in H, 0. We see that the two are 
remarkably alike. Their overall scale agrees because they were fit to the same total energy, 
and they both come from an initial distribution of about the same size, but they were arrived 
at very differently and their agreement is certainly both satisfying and surprising. 



III. COHERENT STATE CALCULATION 

We construct quantum coherent states for the u and tt mesons from the classical fields 
obtained in the previous section using the same methods we used before [0J , . Recall that 
we wish to construct coherent states that have fixed four-momentum and isospin. To do this 
we define a pion field operator that creates pions at the space-time position x and pointing 
in the isospin direction T by 

F(x) = J d 3 kf{k)ai ■ fe~ ik - x , (23) 
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where in the exponent under the integral, the energy component of k is taken as the on-shell 



value k = yk 2 + m 2 . The corresponding omega field operator is given by 

G(x) = J d 3 k'g{k')bi,e- ik '- x . (24) 

with appropriate definition of the fourth component of k! . In these forms, at creates an 

isovector pion with momentum k and bt, creates an isoscalar omega with momentum k' . 

To impose both energy-momentum and isospin conservations, we take the quantum co- 
herent state of the pion and omega system in the radiation region to be given by 

\I,I Z ,K>= / a X a ' e^Y IIz (f)(e F ^ +G ^ - 1 - F(x) - G(x))\0 > . (25) 



As we explain in |TJ , we subtract the "one" as well as the one particle (pion or omega) states 
since they are not permitted by energy-momentum conservation and make the calculation 
numerically unstable. The states defined in (EBft are not normalized, but they are orthogonal. 
We have 

< I, I z , K\f, I' z , K' >= 5\K - K')S lf S Izf l(K) (26) 
The normalization factor is given by 

I{K) = I (2^r ^(r)^(T>^(e^W f - f, ^W - 1 - p n (x)f ■ f - p w (x)) (27) 
where 

Prix)= [ d 3 p\f(p)\ 2 e-^ (28) 



and 



p u {x)= I d*p\g(p)\ 2 e-^ (29) 



The normalization integral is difficult to calculate numerically, even after the subtraction. 
Hence we use the expansion method we developed before 0, generalized to the case of two 
meson types. We thus get 

AK)= E I{K n N ^ ) f ^I). (30) 

The probability for finding N n pions and N u omega is given by 

^^^^w^" (31) 

where 

I{K,N n) N w ) = J 5\K-Y,Pi ~ E?i) Ii d3 Pi\f(Pi)\ 2 II Ail<?(?i)| 2 (32) 

i=l 3=1 i=l j=l 
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and 

J I > N„ and / - JV W is odd 

~~ i Mxl I < N and I — N is even ^ ' 

If either A/^ or A^ equal zero, the corresponding terms will not appear in Eq. (|32|) . 
The mean number of 7r's of isospin type \i in the state is given by 

N ^ = ^l pj/™ ^ /.(^y//-(^)t^ e ^Ar(x)(^W^^W - 1), (34) 
and the mean number of cu's by 

^ = ^/ ( ^4^f ^(r)^(T>^>,(x)(e^^+^W - 1). (35) 

These can again be calculated using the expansion method. It is clear that correlations can 
be calculated as before p|, but we will not do so here. 



IV. RESULTS 

Our calculational plan is now clear. We begin with a spherically symmetric "blob" of 
Skyrmionic matter of size fixed by the 7r-meson mass and amplitude fixed by a total energy 
of 2 GeV, but with baryon number zero. We use Skyrme dynamical equations modified to 
include the omega field to propagate the classical pion and omega fields outward into the 
radiation zone. These classical radiation fields are then used to construct quantum coherent 
states for the 7r-mesons and for the cu-mesons, and the coherent states are projected onto 
states of good isospin and four-momentum. We use these states to find mean meson numbers. 
Note that except for the paramters in the Skyrme model, and these are fixed by nucleon 
physics, there are no free parameters in our picture. 

We first calculate mean meson numbers from the purely classical fields. This is completely 
equivalent to calculating the mean numbers in the coherent state without the isospin and 
four-momentum projections. We find 

N w = 4.4 
N u = 0.86 
^total = 7 -° 

where ^otal ^ s ^ ne total number of pions, coming from both direct pions and from omega 
decay. It is given by ^total = ^ + 3A^. Our case corresponds to a percentage of secondary 
pions of 37%. Note that the mean number of omega mesons is quite small. It is usual to 
say that coherent state methods are best applied in the limit of a large number of quanta, 
and 0.86 is certainly not a large number. Nevertheless we use coherent states methods for 
this case both to remain parallel with the treatment of the pions, and because we have little 
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idea of what else to do. An after-the-fact justification comes from the close agreement with 
data. 

We now turn to a calculation of the number of pions and omega meson from annihilation 
with both isospin and four- momentum projections included as in (|34]) and (j35|). In the 
table we show the number of pions (now separated by charge type) and number of omegas 
for each of the isospin channels that can be reached in nucleon-antinucleon annihilation 
at rest. We also show the total number of pions of each charge type as well at the total 
number of all types. These are calculated combining the direct pions with those from omega 
decay using the decay mode u — > 7r + + 7r° + 7r~. We also give the pion variance for each 
channel. We see, as we expect from its large mass, that the mean number of omega mesons 
is decreased by the imposition of energy and momentum conservation. The mean number 
of pions, the variance, and the fraction by charge type is quite close to what is observed 
||. The percentage of secondary pions is 23% for all three isospin channels. This too is 
roughly what is seen when allowance is made for pions from rho mesons, that we have not 
yet included. It may by argued that since the mean number of pions comes out about right, 
there is no room for pions from rho mesons. That argument is incorrect, because it neglects 
the important effect of energy conservation. Note that including omega mesons has not 
substantially raised the mean number of pions from what we obtain in our earlier, pions 
only, calculation Q, @. Because the total energy is fixed, when we add degrees of freedom, 
the energy is redistributed, less goes into direct pions and more into indirect, but the total 
remains about the same. We expect this will be true when rho mesons are included, but 
that then the number of secondary pions will be close to the experimental value of 40%. 

The pion averages are interesting, but the pion number spectrum is a more stringent 
test of the formalism. We are interested in the probability, P n , of finding exactly n pions in 
annihilation at rest, from states of fixed isospin and four- momentum. This can be calculated 
from the joint probability of finding N m pions, and N u omega mesons, p(N n , using the 
relation 

P n = ]T p(N v ,N u ). (36) 

n=N n +3N w 

The joint probability, p{N. K) N u ), is given by Eq.(^l|). This is the generalization to two meson 
types of the method used for pions only in ||. In Figures 11, 12 and 13, we show P n as 
a function of n for each of the three isospin cases. The solid squares are the results of our 
calculation, the open circles are a normalized gaussian distribution fixed to the same average 
and variance as our calculation. We see, as we saw in our previous work HJ], 0, that with 
four-momentum conservation imposed, our calculation is indistinguishable from a gaussian 
distribution, even though our reaction mechanism is very far from a statistical one. The 
pion number distributions in the three figures look very much like the empirical ones |13 



In our previous, pions only, calculation 0, we found that 1 = states could go only into 
even numbers of pions and 1 = 1 states into odd numbers. No such effect is seen in the 
data. Note that by including the omega meson we have smoothly removed any vestige of 
this even-odd effect. All of this gives futher support to the quantized classical wave picture 
of annihilation. 
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V. CONCLUSION 



We have shown that a classical treatment of nucleon-antinucleon annihilation at rest 
based on Skyrme dynamics extended to include the omega meson and then quantized using 
coherent states gives an excellant account of the pion spectra seen in annihilation. This 
treatment also removes an odd-even artifact of the pions only treatment we presented before 
[1J, 0. The treatment presented here also is more firmly based in the classical dynamics 
and hence more rooted in QCD. In fact our treatment now has no free parameters, and 
yet fits the major trends of the data very well. Its major failing is that we find only 23% 
of annihilations going via vector mesons, while empirically that number is closer to 40%. 
This missing fraction is presumably due to annihilations into rho mesons, which are not yet 
in our picture. We argue that including them will boost the fraction of annihilations into 
vector mesons while enery-momentum conservation will prevent that boost from spoiling 
argeement with pion number distributions. 

Our approach to annihilation is part of a widening family of approaches to strong inter- 
action physics problems that fall squarely in the domain of non-perturbative QCD, but that 
involve fairly large energy release. These problems are very difficult to solve in full QCD. 
The alternative is to solve them first in classical QCD (CCD). There the dynamics is difficult 
but often tractable. In our case it is the Skyrme model modified to include omega mesons. 
Once the classical fields have been obtained, they can be quantized in the radiation zone 
using the method of coherent states and some of the important quantum numbers can even 
be imposed. These quantized states can then be compared with experiment. In our case we 
find good agreement, as far as we have gone, and are encouraged to go further. In particular 
we will study pion correlations, annihilation in flight, the effect of the rho meson and related 
phenomena. It should be noted that one of the benefits of the classical dynamics-coherent 
state approach is that all final channels are treated together. Further afield this direction, 
using classical QCD first and quantizing later, is being applied to disoriented chiral conden- 
sates, Centauro events fT3}| , [jT5| , flip] , |I7| , fT8|j , and much else. It may have other uses 



in high energy heavy ions collisions and in other problems of coherent hadronization. 
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Figure Captions 



Figure 1. Pion field configuration F as a function of r and t. 

Figure 2. uq (in units of omega mass m) field as a function of r and t. 

Figure 3. u r (in units of omega mass m) field as a function of r and t. 

Figure 4. Total energy density multiplied by 4nr 2 (in units of MeV/fm) as a function of r 
and t. 

Figure 5. Pion energy density multiplied by 4nr 2 (in units of MeV/fm) as a function of r 
and t. 

Figure 6. Energy density associated with u r multiplied by 4irr 2 (in units of MeV/fm) as a 
function of r and t. 

Figure 7. Interaction energy density multiplied by 4irr 2 (in units of MeV/fm) as a function 
of r and t. 

Figure 8. Pion momentum distribution. The horizontal axis is the scaled momentum in 
units of inverse omega mass. The vertical axis is the dimensionless reduced pion 
distribution function. The area under the curve gives the total direct pion number N n . 

Figure 9. Omega momentum distribution. The horizontal axis is the scaled momentum in 
units of inverse omega mass. The vertical axis is the dimensionless reduced omega 
distribution function. The area under the curve gives the total omega number N u . 

Figure 10. Pion momentum distribution. Solid line: Skyrme calculation; Dot-Dashed Line: 
Phenomenological pion momentum distribution used in [1,2]. The horizontal axis is the 
scaled momentum in units of inverse omega mass. The vertical axis is the dimensionless 
reduced pion distribution function. The area under the curve gives the total pion number 
n. 

Figure 11. Pion number distribution in / = 0, I z = channel. Solid squares are given by 
the distribution P n . Circles are the gaussian distribution with the same mean and variance. 
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Figure 12. Pion number distribution in / = 1, I z = channel. Solid squares are given by 
the distribution P n . Circles are the gaussian distribution with the same mean and variance. 

Figure 13. Pion number distribution in / = 1,I Z = 1 channel. Solid squares are given by 
the distribution P n . Circles are the gaussian distribution with the same mean and variance. 
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TABLES 



TABLE I. pion number distribution. N n is the total direct pion number and is the total 
omega number, no is the total 7r° number including those from omega decay. Similarly n + and n_ 
are the total ir + and tt~ numbers including those from omega decay, n is the total pion number 
and a is the standard deviation for n calculated from the distribution function P n . 



Channel 


N n 




no 


n + 


n_ 


n 


a 


7 = I z = 


4.94 


0.60 


2.25 


2.25 


2.25 


6.75 


0.79 


7 = 1 7 2 = 


5.17 


0.48 


3.99 


1.31 


1.31 


6.61 


0.84 


7 = 1 I z = l 


5.20 


0.48 


1.32 


3.16 


2.16 


6.64 


0.84 



15 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405351vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405351vl 



This figure "fig3-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405351vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405351vl 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405351vl 



This figure "fig3-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405351vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405351vl 



This figure "fig2-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405351vl 



This figure "fig3-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405351vl 



This figure "fig3-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405351vl 



This figure "fig3-5.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405351vl 



This figure "fig3-6.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405351vl 



This figure "fig3-7.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405351vl 



